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ABSTRACT 

Using high resolution SPH simulations in a fully cosmological ACDM context we study 
the formation of a bright disk dominated galaxy that originates from a "wet" major merger at 
z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early 
major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly 
following the last major merger and the present day properties of the final remnant are typical 
of early type spiral galaxies, with an i band B/D ^ 0.65, a disk scale length of 7.2 kpc, g — r 
= 0.5 mag, an HI line width (W20/2) of 238 km/sec and total magnitude i = -22.4. The key 
ingredients for the formation of a dominant stellar disk component after a major merger are: i) 
substantial and rapid accretion of gas through cold flows followed at late times by cooling of 
gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation 
during mergers and iii) relative fading of the spheroidal component. The gas fraction of the 
progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of 
gas from the local environment plays a major role in the regrowth of disks and in keeping 
the galaxies blue. The results of this simulation alleviate the problem posed for the existence 
of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at 
relatively low z. 

Key words: galaxies: formation, evolution, interactions, methods: N-Body Simulations. 



1 INTRODUCTION 

Within the ACDM framework the build up of galaxies and 
their parent dark matter (DM) halos occurs t hrough a series of 
mergers and accretion of more diffuse matter jprenk et"ai]|l985l : 
[Springel & et. al. 2005). The classic models of galaxy formation 
(White &Rees 1978; Fall 1983; White & Fren k 1991.) an d sub- 
sequent work iD alcanto n et all 119971; [Mo et al] [1993; ISilkll200ll ; 
iBower et all200Q ; iZheng & et al.l2007l ; ISomerville et alj2008h as- 



sume that gas infalls inside the parent dark matter halo and subse- 
quently cools to temperatures low enough to fragment and form 
stars. The morphology of galaxies and the fraction of stars in 
their disk and spheroidal components are set by a combination 
of competing physical processes, as numerical simulations have 
shown that violent relaxation in mergers between similar size 
galaxies can destroy or dyna mically heat their ste llar disks and 
turn them into spheroids l iBarnes & HemguisJl 19961) . On the other 
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hand, remaining and subsequently accreted gas can regrow a disk 
teaugh et al] 1 19961 ; Isteinmetz & NavMrol l2002h if given enough 
time. It is then a prediction of hierarchical models that the morphol- 
ogy of a galaxy is not assigned "ab initio" but it might change of- 
ten over cosmic times, as the spheroidal to disk light ratio changes 
back and forth. Observed merger remnants have disky isophotes 
l|Rothberg & Iosepri'2004') or evidence for young disk components 
(McDermid & et al. 2006.), hinting at disk regrowth and supporting 
the above scenario. 



Observationally, major mergers are observed to be common in 
the redshift range - 1 . Their number density ha s been estimated to 
be of the order of 1 0~<^~'^>h~^Mpc~'^ Gyr~^ , jLin & et alj|20o3; 
llogee & et al.ll2008l) . possibly increasing with redshift. As numer- 
ical work suggests that a consistent supply of gas is important to 
regrow a disk component ( Robertson et al. 2006), it is relevant that 
observational evidence for mergers between actively star forming 
galaxies has also emerged, showing that at redshift ~1, 70% of the 
merging pairs are between blue (or "wet", defined as rest frame 
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g-r <0.65) pairs; however the fraction of blue and likely gas rich 
mergers decreases by the present time jLin & et alJlOOSi) . 

Numerical work has also highlighted a possible problem 
for the existence of disk galaxies at low redshift, showing 
that interactions, mergers and accretion events are common at 
every redshift for galaxy sized h alos in ACDM models 
jMaller et all 12004 IStewart et al.' '20081) and potentially destruc- 
tive for disks jlbth & Ostriker 1 99a iKazantzidis et alj |2007| : 
iBullock et al]|2008l ; rPurcell et al.ll2d08l) . This makes it potentially 
difficult to reconcile the observed present day population of disk 
dominated galaxies with ACDM, if a quiet merging history is 
essential for their existence. But do disk galaxies really require 
a quiet merging history to be observed as such at the present 
time? What are the time scales of disk growth and destruction? 
How quickly do disk galaxies regrow after a merger? Theoret- 
ical models also highlight the difficulties in predicting the out- 
comes of galaxy mergers and the subsequent regrowth of stellar 
disks. These difficulties arise from several factors. Given the de- 
coupling of the baryonic cores from the parent halos it is diffi- 
cult to predict robust merging rates of galaxies (as compared to 
their DM halos) unless cosmologic al simulations with hydrody- 
namics are used ( Mailer et al. 2006*). Also, numerical simulations 
have shown how the bulge to disk ratio (B/D) resulting from a ma- 
jor merger might depend on the orbital parameters, the internal spin 
of the progenitors, the gas fraction of the p arent dis ks and the effi- 
ciency of SN feedback feames & Herngui st 1996; C ox et alj2003 ; 
IScannapieco et alj2008 ). 

Recent observational and theoretical work has pointed out 
ways to alleviate the possible problem of a h igh m erger rate for 
the survival of disk galaxies. iHopkins et all ( l2008h showed that 
angular momentum loss of the gas component is not necessarily 
catastrophic even in 1:1 mergers. Disks can survive or rapidly re- 
grow, provide d that the gas fraction in the disks of the progen- 
itors is high ([Robertson et alj |2006| ; [Robertson & BuUockl |2008| ; 
iBuUock et al.ll2008h . If f eedback is a ble to suppress star formation 
during the merger event |Brook et ah 2004 ; Robertson et al. 2006; 
IZavala et alj2008l ; lGovemato et al.ll26o7l) the existing cold gas can 
settle on a new disk plane and start regrowing a stellar disk. 

Hydrodynamical simulations indeed show that cold flows 
(cold gas that flows rapidly to t he center of galaxies from fila- 
menta ry structures around halos jKeres et"aill2005l ; lDekel & et al.l 
1200^)) 1 play a major role in the build up of disks in galaxies 
dBrooks & et. alluOOa ). Gas accreted through cold flows arrives to 
the central stellar disk on a time scale a few Gyrs shorter than 
gas that is first shocked to the virial temperature of the host halo 
and then cools onto the disk, leading not only to early disk star 
formation, but the creation of a large reservoir of cold gas. No 
matter its origin, late infalling gas would likely have a higher an- 
gular momentum conte nt than material accreted at earlier times 
dUuinn &Binnevl[l992l) and gas in the merging disks would ac- 
quire a coherent spin set by the orbital parameters of the binary sys- 
tem and the internal spins of the parent galaxies. A feedback+cold 
flows model is particularly attractive as the mechanism for the sur- 
vival/regrowth of gas (and then stellar) disks as it provides a natural 
explanation to the fact that more massive galaxies tend to ha ve large 
B/D ratios teenson et all l200l iGraham & WorlevI |2008|) . while 
smaller galaxies are likely disk dominated. At large galaxy masses 
energy feedback from supemovae (SNe) becomes ineffective at 
suppressing star formation while cold flows become inefficient at 
carrying a supply of fres h gas necessary to regrow a stellar disk 
dDekel & Bimboin]|2006h . Combined with th e relative higher fre - 
quency of major mergers for massive galaxies touo & Whitel2008l) 



this framework leads to the build up of a larger stellar spheroid and 
disfavors the quick rebuilding of stellar disks in massive systems. 

Only recently have numerical simulations been able to di- 
rectly co mpare the observable quantitie s of t he outputs with 
real data djonsson 2006; Chakrabarti et a Li i2007l ; iLotz et al.ir2008l; 
ICovington et alj|2008l ; feocha et alJi2008l) . rather than simply pre- 
dicting the mass distribution of the simulated stellar systems. This 
is a crucial point as the stellar spheroids will fade drastically after 
a major merger, while reforming disks will contain a large fraction 
of younger and brighter stars. Different mass to light (M/L) ratios 
for the two components will skew the observed photometric light 
ratios compared to the underlying stellar masses. 

Despite all of this progress, the effect of mergers on the ex- 
isting disk components, and of feedback and cold flows on the re- 
growth of stellar disks have not been studied in detail in fully cos- 
mological simulations of major mergers. Here we present results 
from a fully cosmological, smooth particle hydrodynamic (SPH) 
simulation where late gas accretion plays a major role in the rapid 
regrowth of a dominant stellar disk in an L* galaxy after a z =0.8 
major merger. In this work we focus on the physical processes that 
drive the regrowth of the disk as identified by its kinematical prop- 
erties, but we also measure the structural properties of the simu- 
lated galaxy based on the light distribution, i.e. in a way closely 
comparable to observations. This result helps solve the apparent 
contradiction that strong interactions at relatively low redshifts are 
common even for DM halos that are likely to host bright disk galax- 
ies. The paper is organized as follow: §2 discusses the simulations, 
§3 describes the evolution of the system, §4 describes the observa- 
tional properties of the merger remnant and §5 the assembly of the 
disk component. The results are then discussed in §6. 



2 DESCRIPTION OF THE SIMULATION 

The simulation described in this paper is part of a campaign of 
high resolution simulations aimed at studying the formation of 
field galaxies in a WMAP3 cosmology (Slo=0.24, A=0.76, h=0.73, 
(78=0.77, S7i,=0.042). Our sample of halos has been selected from 
low resolution volumes of size 50 and 25 Mpc (the latter for dwarf 
size halos) using the "zoom-in" tecnique to ensure a proper treat- 
ment of tidal torques iKatz & White|[l993l) . The sample covers two 
magnitudes in total mass (from 2x 10^° to 2x lO'^'^ Mq) sampling a 
representative range in halo spins and epochs of last major merger. 
The galaxy described in this paper has a total mas s of 7x1O"M0, 
spm A = 0.04 (as defined in lBullock et"ai] ( 1200 ih ) and a last ma- 
jor merger at z=0.8. The environment density at z=0 is typical of 
field halos with an over density 5 pip = 0. 1 (0.2) measured over a 
sphere of radius 3 (8) h~^ Mpc. The final galaxy has 650k, 320k 
and 1.8m DM, gas and star particles within the virial radius of 
the galaxy at the present time, with particle masses: 1.01x10^, 
2. 13 X 10^ and 6.4 x lO"* Mq for each DM, gas and star particle at 
the moment of formation. The force spline softening is 0.3 kpc. 
The minimum smoothing length for gas particles is 0. 1 times the 
force softening. The si mulations were perfo rmed with the N-body 
SPH code GASOLINE dWadslev et al.l2004h with a force accuracy 
criterion of 6* = 0.725, a time step accuracy of ?7 = 0.195 and a 
Courant condition of rjc = 0.4. The adopt ed star format i on an d SN 
schemes have been described in detail in Stinso n et ai] d2006t) and 
Governato et al.l (2007). Our "blastwave" feedback scheme is im- 
plemented by releasing energy from SN into gas surrounding young 
star particles. The affected gas has its cooling shut off for a time 
scale associated with the Sedov solution of the blastwave equation. 
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Figure 1. SDSS rest frame, unreddened i-band surface brightness images of the merging system at three different stages. Left (panel a): z=1.0 close 
to the first pericentric passage. Center (panel b): the merger remnant at z=0.6, IGyr after the stellar cores of the two galaxies have coalesced. Right (panel c): 
the system observed at the present day: the halo has faded considerably while the galaxy has grown an extended thin disk. The i-band disk scale length R^j is 
7.2 kpc and the light ratio of the kinematically identified Bulge and Disk components (B/D) is 0.49, while a GALFIT 2D decomposition gives a B/D = 0.65. 



which is set by the local density and temperature of the gas and 
the amount of energy involved. At the resolution of this study this 
translates into regions of ^ 0.2-0.4 kpc in radius being heated by 
SN feedback and having their cooling shut off for 10-30 million 
years. The effect is to regulate star formation in the disks of massive 
galaxies and to lower the star formation efficiency in galaxies with 
circul ar velocity in the 50 < Vc <150 km/sec range jSrooks et al.l 
l2007h . At even smaller halo masses (Vc < 20-40 km/sec, with Vc 
= sqrt(G X M/r)) the collapse of baryons is partially suppresse d 
by the cosmic UV field ( Hoeft et al. 20od ; lGovernato et alj|2007h . 
here modeled following iHaardt & Madai] The simulation 

applied a correction to the UV flux fo r self shielding of dense gas 
as introduced in Pontzen et al .1 Igoosh and low temperature metal 
cooling jMashchCTko et a l. 2008). It is important to note that the 
only two free parameters in the SN feedback scheme (the star for- 
mation efficiency and the fraction of SN energy coupled to the ISM) 
have been fixed to reproduce the properties of present day galax- 
ies (star formation rates, Schmidt law, cold gas turbule nce, disk 
thickness) over a range of masses dGovernato et al.ll2007l) . Without 
further adjustments this scheme has been prov en to reproduce the 
relation between metal licity and stellar mass terooks et allbOOTl : 
iMaiolino & et al.ll2008h and the abu ndance of Damped Lyman a 
(DLA) systems at z=3 jPontzen et al . 2008). However, even at this 
resolution the central regions (r<l-2 kpc) of galaxies remain still 
partiall y unresolved and likely form bulges that are too concen- 
trated dGovernato et al. I l2008l : lMaveretalJl2008l) . In this respect 
the mass of the bulge component of our simulated galaxy has to 
be considered an upper limit imposed by current resolution limits, 
especially at high z, where the number of resolution elements per 
galaxy is less. 

To properly compare the outputs from the simulation to real 
galaxies and make accurate estimates of the observable properties 
of gal axies, we used the Monte Carlo radiation transfer code 56W- 
RISE ('Jo nssonll2006l) to generate artificial optical images (see Fig- 
ure 1) and spectral energy distributions (SEDs) of the outputs of our 
run. SUNRISE allows us to measure the dust reprocessed SED of 
every resolution element of our simulated galaxies, from the far UV 
to the far IR, with a fully 3D treatment of radiative transfer. Filters 



Table 1. Merger remnant propeilies at z=0.8 (shortly after the merger) and 
at the present time. is the disk scale length and B/D is the bulge to 
disk light ratio measured for the kinematically identified components. (B/D 
ratios in parentheses are measured using a 2D photometric decomposition 
with GALFIT of the face on projection of the light distribution). Total mag- 
nitudes and colors have been measured in SDSS filters, including the effects 
of dust (disk inclination 45deg) and in the AB system. 





z=0. 


z=0.8 


Tot Mag i band 


-22.3 


-22.7 


R^ i band (kpc) 


7.2 


4 


B/D i band 


0.49 (0.65) 


1.1(1.4) 


B/D (stellar mass) 


0.87 


1.16 


g -r 


0.5 


0.4 


SFR 


2.2 Mg/yr 


6 MQ/yr 


disk stellar mass 


3.24x1010 


2.07x1010 



mimi cking those of the SDSS survey dAdelman-McCarthv & et aT] 
l2006l) are used to create mock observations. 

To measure the rotational velocity of the remnant and other 
galaxies in the sample in a way comparable with observations we 
used the spatial and kinematic distribution of cold gas in the disk of 
our simulated galaxy. Using the HI fraction calculated by GASO- 
LINE, we determined the HI line width, (W20), by finding the 
width of the HI velocity distribution at 20% of the peak. The value 
of W20/2 is then used as a measure of the galaxy rotation velocity 
at different redshifts. This measurements reflects the mass weighted 
position-velocity distribution of cold gas inside the central region 
of the galaxy and in observed bright galaxies it is usually asso- 
ciated with the peak rotational velocity. Because the mass inside 
the c entral few kpc is li kely overestimated owing to resolution ef- 
fects dMaver et alj2008l) , VK20/2 provides a slightly larger measure- 
ment of the rotational velocity than that obtained from the rotational 
speed of young stars at 2.2 or 3.5 disk scale lengths (typically 10- 
20kpc for bright galaxies), and thus is a useful upper limit for a 
comparison with real data. 
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Figure 2. SFR vs time for the last major merger progenitors (red dotted and 
blue dashed) and for the whole galaxy (solid) . The peaks of the dotted and 
dashed curves correspond to the two major mergers between the original 
four progenitors. 
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Figure 3. Lower Panel: Time evolution of the cold disk gas (T<4x 10* K) 
scale length and the stellar disk scale length (measured in the unreddened i 
band) of the final merger remnant. Upper Panel: The B/D ratio for the final 
remnant as a function of redshift: dashed: stellar masses. Dotted: B/D light 
ratio (unreddened i band) with a GALFIT decomposition based on the 2D, 
face on light profile. Continuous: B/D light ratio (unreddened i band) of the 
kinematically defined disk and bulge components. 



3 DYNAMICAL AND PHOTOMETRIC EVOLUTION OF 
THE MERGING SYSTEM 

The galaxy studied here was singled out for its particularly inter- 
esting assembly history, as the build up of its stellar component in- 
volves numerous major mergers, seemingly a hostile environment 
to build a significant stellar disk. At z=3 the four most massive 
progenitors form a hierarchy of two binary systems roughly aligned 
on the same large scale filamentary structure. Each halo has a to- 
tal mass ~ 7x 10^'^ Mq and has formed a rotationally supported 
stella r disk fed by stro ng cold flows, typical of galaxies at that red- 
shift te rooks & et. allboOS.) . In this work we define as "cold flow" 
gas that has never been shocked to 3/8 the virial temp erature of the 
parent halo jKeres et alJl2005l . lBrooks & et. allllOOa) . The average 
cold (T<4xl0'' K) gas fraction in the disks of the four galaxies 
is 25% (defined as the fraction of total baryons in the disk). Both 
pairs merge by z=2. Just before the mergers, the four progenitors 
have rest frame B magnitudes in the range -21.1 — 20.2 (-2L5 - 
-20.6 in the r band) and rest frame g ~ r colors around 0.3 - 0.4 
(unless specified all global magnitudes and colors in this work are 
in the rest frame AB system and include the effects of dust red- 
dening measured at a 45 deg inclination). These two early mergers 
are then "wet" i.e. between galaxies with blue colors as defined in 
lLin&etal.U2008h . 

Both merger remnants quickly reform extended gas disks from 
a combination of freshly accreted gas and gas already in the progen- 
itors disks that was not turned into stars during the merger. For each 
merger the star formation history (SFH) peaked at 18 and 32 Me/yr 
respectively (Figure 2). At z=1.6 (2 Gyrs before the final merger) 
the two disk galaxies formed from the early mergers have again 
very similar magnitudes: -2L6 (-22.1) and -21.3 ( -21.8) in the B (r) 
band respectively, close to the B band L* at that z jMarchesini et al.l 
|2007). Their p — r color = 0.4, makes the m bluer than most galaxies 
of similar brightness at the present time jLin & et al.ll2008l) . 

The final major merger of these two progenitor galaxies be- 
gins at around z=l when their dark matter halos, flowing along 
the same filamentary structure, first overlap. The galaxies plunge 
in on fairly radial orbits, with the internal spins of the two disks 
roughly aligned with the orbital angular momentum vector (Fig- 
ure la). Note that this is not necessarily a configuration favorable 
to the survival of gaseous disks, as noted by Hopkins et al. (200^. 
After two close passages, the two galaxies coalesce by ~ z=0.8, 
i.e., 1 Gyr after the merger commenced. The mass ratio of the 
merging halos is 1.2:1 while the i band brightness ratio is 1.6:1. 
During the merger, the global star formation rate of the system 
is enhanced and peaks at 11 Mq with subsequent star formation 
rates in the remnant dropping rapidly by almost factor of three 
to ~ 4 M0/yr (Figure 2). The moderate star formation (SF) en- 
hancement is consistent with estimates for i nteracting systems in 
the same redshift range ( I Jo gee & et al.ll200^ . Given the galaxies' 
properties outlined above, even this final merger is then clearly 
identified as "wet". However, the disks of the two galaxies have 
gas fractions around 20%, so they are only relatively gas-rich com- 
pared to the present day popul ation of galaxies of similar brightness 
( Gar cia- Appadoo et al.ll2008|) and have equal or lower gas abun- 
dance as z ~ 2 galaxies terb et alj20M) . 

Once again, following the final major merger a gas disk 
rapidly regrows and star formation is mainly concentrated in the 
reforming stellar disk. The galaxy remains relatively unperturbed 
after z=0.8. Shortly after the merger the spheroidal component of 
the newly formed galaxy dominates the light distribution (Figure 
lb), although a disky component is already visible. The g — r color 
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Figure 4. The time evolution of the angular momentum content per unit 
mass of the kinematically identified z=0 baryonic disk (gas and stars) and 
the dark matter of the merger remnant. The reference frame is defined by 
the center of mass of the selected particles. Disk material conserved most of 
the angular momentum gained by z=1.5. The vertical line marks the epoch 
when the halos and the baryonic cores merge during the last major merger. 

is 0.5 and remains stable to the present time. We verified that the 
disky component at z=0.8 is indeed associated with a thick stellar 
disk supported by rotation. By z=0 the galaxy has regrown an ex- 
tended thin disk, and while about 50% of the baryons in the halo 
have been turned into stars, the spheroid has faded considerably. 
The halo has faded in the reddened B band by 0.8 mag to B = - 
21, due to the aging of the spheroidal component. The disk clearly 
dominates the light distribution in the SDSS i and bluer bands (Fig- 
ure Ic). 

These results show that at all stages the progenitors of the fi- 
nal galaxy can be identified with a normal population of moderately 
gas rich disk galaxies with colors and SF rates comparable with the 
galaxy population observed in the 3>z>l range. It is also impor- 
tant to emphasize that the cold gas fraction in all the disk progeni- 
tors suggests that idealized initial conditions with disk gas fraction 
as high as 50% are not a necessary condition to have "wet" mergers 
and rapid regrowth of stellar disks. As discussed in §5 the funda- 
mental ingredient for disk regrowth and blue colors is the fast re- 
fueling of gas from the hot halos and the surrounding cosmic web, 
highlighting the necessity of a fully cosmological approach to the 
problem of the abundance of galaxy disks^. 



4 PROPERTIES OF THE MERGER REMNANT FROM 
Z=0.8 TO THE PRESENT 

In this section we decompose our simulated galaxy into its kine- 
matic components and "observe" some of its properties at the 
present day (they are summarized and compared with those of the 
galaxy at z=0.8 in Table 1). Does the galaxy have the photometric 
and kinematic properties of disk dominated galaxies? 

^ Movie at www.astro.washington.edu/fabio/movies/Merger.mpg 
© 2002 RAS, MNRAS OOO.fTlflOl 



At the present time the disk clearly dominates the light distri- 
bution even in the relatively red i band (Figure Ic) The galaxy has a 
total magnitude is i —-223 (-22.4 unreddened) mag and global red- 
dened color g — r =0.5, consistent with those of luminous present 
day disk galaxies l lLin & et alj|2008l) . 

To measure its morphology in a quantitative way the different 
components of the galaxy were first identified using their kinematic 
and spatial information and classified as bulge, halo and disk. This 
is a crucial step to relate each galaxy component to its physical ori- 
gin. First the disk plane is defined using the cold gas in the central 
few kpc of the galaxy, then disk stars are defined as stars whose 
specific angular momentum perpendicular to the disk plane (jz) is 
a significant fraction of the maximum angular momentum of a cir- 
cular orbit with the same binding energy (jc), i.e j^/jc > 0.8. Parti- 
cles on circular but inclined orbits (more than 30 deg) are excluded. 
Bulge and halo stars were then identified based on their radial or- 
bits and their binding energy (bulge stars being more bound). The 
energy separation criteria between halo and bulge stars corresponds 
to the radius at which the spheroid mass profile changes slope (halo 
stars having a shallower profile than bulge stars) and in our simu- 
lated galaxy sample separates an older and metal poor population 
(the halo) from bulge stars that are more metal rich. Halo stars con- 
tribute a fraction (~I5%) of the total stellar mass within the virial 
radius of the galaxy, but the halo central density is two orders of 
magnitude lower than that of the bulge. Hence the details of the 
bulge/halo decomposition do not change our conclusions. At z=0 
the kinematically identified disk, bulge and halo stellar masses are 
respectively: 3.4, 2.7 and 1.x IO^^Mq. We then imaged each sep- 
arate component using SUNRISE and measured their structural pa- 
rameters using the unreddened images. We focused on a structural 
analysis of the unreddened components, avoiding the additional 
layer of complexity given by the details of the dust distribution, 
which will be explored in future papers with a larger number of 
galaxies. However, we have verified that our findings do not change 
if the reddened images are used instead. 

How and when did the disk reform after the last major merger 
at z=0.8? The stellar disk and bulge components were identified at 
different redshifts after the last major merger event. To better under- 
stand the B/D ratio evolution of the merger remnant we measured 
B/D in three different ways. We used the kinematic decomposition 
to find a) the stellar mass ratio of the bulge and disk components 
and b) their relative flux ratio in the i unreddened band. Then we 
analyzed the unreddened, face-on 2D light distribution created by 
SUNRISE using all th e galaxy star p articles (including halo parti- 
cles) with GALFIT (Peng et al. 2002) to find c) the B/D ratio as 
determined by a fit to the surface brightness profile. Figure 3 shows 
how the i band disk scale length and the B/D ratio evolve with time. 
Shortly after the merger event the disk component is already visi- 
ble edge-on, then the bulge component fades relative to the disk and 
the disk becomes more extended. At z<0.4, or about 3.5 Gyrs after 
the merger, the disk dominates both in terms of the light contribu- 
tion and the B/D ratio decreases further by z=0. All measurements 
agree on the same trend of B/D decreasing with time. At the present 
time the disk extends almost far as 20 kpc in radius from the galaxy 
center and the stellar disk scale length Rd is 7.8 and 7.2 kpc in the 
B and i bands, respectively (Figure 3), consistent with observations 
of real galaxies that show larger Ra in bluer bands. Smaller B/D 
ratios are obtained using the light distribution, more sensitive to the 
younger ages of disk stars. At z=0 the B/D stellar mass ratio of 
the kinematically defined components is 0.87, but the unreddened i 
band light ratio is only 0.49. GALFIT shows the steepest trend with 
age, and shortly after the merger it underestimates the disk compo- 
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nent, if by less than 20%. GALFIT gives a fairly precise estimate 
of the light weighted B/D ratio when the stellar disk becomes dom- 
inant. 

Soon after the merger the cold gas disk increases its size by 
nearly a factor of two as the new infalling gas settles on high an- 
gular momentum orbits. The angular momentum of the present day 
disk baryons is mostl y acquired at high z (F igure 4) as predicted 
in analytical models dOuinn & Binne 1 19921) with a fraction of it 
transferred to the DM halo during the last major merger. As this 
gas is gradually converted into stars the stellar disk also grows in 
size, however the cold gas disk remains significantly more extended 
than the stellar one. Shortly after the merger the (unreddened i 
band) disk scale length is 4 kpc, at the present time it is almost 
twice as large, with only minor warping (Figure 7). Most likely 
due to the late assembly of its younger component, this is a fairly 
extended disk for galaxies of this mass, more extended than many 
disks formed in cosmological simulations where the assembly his- 
tory of the galaxy was more "quiet". Encouragingly, we verified 
that the structural properties of the bulge and disk do not change 
much if they are measured using GALFIT on the global unred- 
dened 2D light distribution, i.e. without a prior knowledge of the 
kinematic decomposition. With GALFIT the i band Rd is 6.8 kpc, a 
5% difference. GALFIT finds systematically larger B/D ratios, but 
only by 20% or less. However, B/D ratios decrease if the GALFIT 
fitting is done on dust reddened images. These B/D ratio and are 
quite typical of bright Sa and Sb galaxies that typically have dust 
corrected B/D ~ 0.5 and red Rd ^ 3-7 k pc ( IGraham & WorlevI 
|2008|; [Driver et aI.ll2007ljBenson et al.]|2007h . 

To show how the simulated galaxy of this study relates to the 
general population of real disk galaxies on the Tully Fisher (TF) re- 
lation, we compared its HI velocity width K^2o/2 (238 km/sec) and 
total rest-frame unreddened i-band relation with the disk dominated 
galaxies in our s imulated sample an d with the sample of disk galax- 
ies described in iGeha et al] ( l2006h and references therein (Figure 
5). We find that the agreement between the observed TF and our 
set of simulations is quite good, as the combination of high res- 
olution and the feedback adopted in our simulations yields a good 
match to real galaxies over a wide range of magnitudes and circular 
velocities. The merger remnant object of this study lies well within 
the observed scatter of both the observed and simulated samples, 
confirming that its structural properties are similar to those of typ- 
ical bright spiral galaxies. We will present the scaling properties 
of the full data-set of simulated galaxies in a forthcoming paper. 
An interesting feature of the Tully Fisher plot is the evolution of 
the remnant on the TF plane: the very limited growth of the bulge 
stellar mass (only a few %) and the inside out growth of the disk en- 
sures that the amount of mass within the central region of the galaxy 
does not change, hence W20/2 does not evolve strongly, while the 
overall fading of the stellar components makes the galaxy dimmer 
by less than half a magnitude between redshift 0.8 and (see Ta- 
ble I). Even if the assembly history of our galaxy is not typical of 
galaxies of similar total mass, this result is consistent with observa- 
tions that find a small evolution in the observed galaxy TF relation 
up to z= 1 ( Conselice et al. 2005) along with strong size evolution 
jTruiillo et al.l2007l) . This result is not trivial and we plan to extend 
this analysis to our full sample of simulations in a future work. 

We also verified that the metallicity of the cold gas is 
consistent with the observed stellar mass - metallicity relation 
(8.5 < 12-l-log(0/H) < 8.9, depending on the aperture used, 
iTremonti & et al.ll2004l ; iBrooks et al.ll2007h . An average metallic- 
ity consistent with real galaxies is an important test of the realism 
of this simulation and makes the estimates of the galaxy colors of 
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Figure 5. The location of the simulated galaxy in this work on the present 
day "Tully Fisher" relation, i.e., the unreddened i ban d vs W20/2. Gre y 
dots: Galaxies from a compilation of obsei-vational data iGeha et ai]|2006l) . 
Triangles: other disk galaxies from our simulated sample at similar resolu- 
tion. Filled dot in main panel: the galaxy of this study at z=0. Filled dots in 
upper left insert: different snapshots of the merger remnant from z=0.8 to 
the present. Circular velocities are the W2Q/2 velocity widths from the HI 
distribution of the simulated galaxies. The insert shows the weak redshift 
evolution of the merger remnant on the TF plot. 



the final remnant and its progenitors more robust. Finally, the satel- 
lite system of the remnant includes 1 1 resolved luminous satellites 
within the virial radius of 230 kpc. The faintest has AB B mag = 
-8.7, the brightest -17.9. By tracking the satellites through different 
outputs after the merger event we verified that the galaxy disk un- 
dergoes several fly-bys by small dark satellites, but no significant 
accretion of luminous satellites after the final merger. Many faint 
satellites have undergone severe tidal stripping of both their DM 
halos and of their stellar component. This analysis quantifies the 
dramatic regrowth of the disk component and how, coupled with 
the fading of the bulge and halo components it leads to the forma- 
tion of a galaxy dominated by an extended disk. It also highlights 
the difference between evaluating a galaxy morphological type us- 
ing the mass distribution compared to the light distribution. It is 
encouraging however, that similar trends (growing disk size and 
decreasing B/D ratios) are recovered using complementary tech- 
niques, as GALFIT provides a decomposition into bulge and disk 
components quite similar to that obtained using the full spatial and 
kinematic information. These results support the notion that the ob- 
served B/D ratios in bright galaxies are indeed representative of the 
underlying dynamical disk and spheroidal components. 



5 THE ASSEMBLY OF THE GALAXY DISK 

A detailed analysis of the gas accretion history was performed by 
tracking backward every gas particle that was ever within the virial 
radius of the simulated galaxy and its progenitors. Every star par- 



© 2002 RAS, MNRAS 000,[TlQo] 



Formation of a large disk galaxy 7 



•■*VT"l''J""-" !■ Ji ■ ■ 





Figure 6. Top Panel: The total gas accretion history of the merger remnant, 
divided by the thermodynamical history of the gas. (at high z the lines are 
the sum of the contributions from individual progenitors). Bottom Panel: 
The SFH of disk star particles separated by the history of their parent gas 
particles. The disk stars are kinematically identified at z=0. Dashed: from 
clumpy gas accretion, dot dashed: from unshocked gas (or "cold flows"), 
solid: shocked gas. A significant fraction of stars formed at high z and were 
identified as part of the disk even after two major mergers. At the present 
time they form the thick disk component (see also Figure 7). Star formation 
in the disk is partially disrupted during the last major merger, and stars 
formed from gas cooled from the hot halo form a significant fraction of the 
younger disk component. 

tide comes from a gas progenitor which is uniquely identified at 
the moment of the star particle formation. The thermodynamical 
history of each gas particle progenitor was then studied. This tem- 
perature history was done for each of the original four progeni- 
tor galaxies, then again for the resulting two binary galaxies, and 
finally for the merger remnant. Gas particles were then classified 
as a) "cold flows" if it never belonged to a progenitor halo before 
being accreted onto one of the progenitor galaxies or the merger 
remnant and it never shocked to 3/8 the virial temperature of the 
galaxy, b) "shocked" if the gas was first shocked to > 3/8 the virial 
temperature of the main halo and then cooled onto the disk, and c) 
"clumpy" if the gas particle was assigned to a subhalo before be- 
coming part of the disk (i.e., a halo other than one of the four main 
progenitors and their subsequent merger remnants). 

Figure 6 shows (top panel) the total gas accretion rates to pro- 
genitor galaxies and merger remnant, and (bottom panel) the SFH 
of the resulting disk (as identified at z=Q based on the stellar kine- 
matics), separated by the different histories of each gas particle pro- 
genitor. It is remarkable that a large fraction of the disk formed at 
redshift one or earlier, when the galaxy had not undergone its last 
major merger yet. This result supports the suggestion that major 
mergers do not alwa ys completely destroy the pre-existing disks 
jHopkins et alJizOOa) . although they heat them quite substantially, 
making them thicker. Cold flows build the largest fraction of the 
disk up to low redshift. This analysis shows how the build up of 
the disk of bright galaxies differs from the simple picture of gas 
cooling from the hot halo and assembling into a disk; only at later 
times does gas cooling from the hot halo phase plays a role. From 
Figure 6 it is clear that cold flows dominate the early gas accretion 
history of the progenitor galaxies, and is the dominant source of SF 
in the galaxy disk. Because SN feedback regulates SF in the disk, 
a cold gas reservoir is able to develop that sustains SF until the 
present time. The SF rate drops significantly when the two progen- 
itor disk galaxies merge at z^\. After this time, gas cooling from 



Figure 7. The i band (unreddened) edge on image of the galaxy disk at 
z=0. Upper Panel: fight distribution from stars formed after the last major 
merger. Lower panel: the light distribution from stars formed before the last 
major merger, but still identified as disk based on their angular momentum. 
The image is 40 kpc across The relative brightness of the two components 
is not to scale. 



the hot halo becomes more important, resulting in 30% of the disk 
stars formed in the last 2 Gyrs formed from gas labeled as shocked. 
However, only 10% of the total amount of stars in the disk at the 
present time have formed from gas that was previously heated to 
the virial temperature of the main progenitor and almost 40% of 
the disk stars after z=0.8 formed from cold flows accreted after the 
last major merger. We find that thi s result is quite typica l for disk 
galaxies of total mass < 10^^ M0 terooks & et. afcOOSi) stressing 
the role of cold flows in building early disks and of "shocked" gas 
in rebuilding some of the young, and hence bluer, more metal rich 
and brighter part of the disk. These effects will have to be carefully 
modeled in semi analytical models of galaxy formation. Figure 1 
shows that as the disk grows in mass, most of the stars are added 
to a thin and dynamically cold part of it. To quantitatively evaluate 
this process we separated the stars previously dynamically identi- 
fied as disk into two populations: those that formed before the z= 
0.8 merger and those formed after, i.e in the last ~ 6 Gyrs. Figure 
7 shows the edge-on i band surface brightness image (unreddened) 
of the two populations. Their properties are strikingly different: the 
younger stars form a thin exponential disk with i band scale length 
of 8.4 kpc (measured using stars in the 2-13 kpc range) and scale 
height of just 0.5 kpc (measured as a Gaussian distribution at a ra- 
dius of 8.5 kpc). This scale height is likely as small as allowed by 
the force softening of 0.3 kpc. The older component, while being 
comparable in mass, is 1.3 mag fainter in the i band and 2.2 mag 
in the B band. It also has a much shorter disk scale length, 3.1 kpc, 
and is considerably thicker, 2.7 kpc, again at an 8.5 kpc radius. The 
flaring in the old disk is due to the kinematic selection criteria, that 
exclude particles on circular orbits that are highly inclined from the 
disk plane. 

While beyond the scope of this paper, we speculate that the 
mix of the two components could easily be identified as a thick 
and thin di sk components if observed using typical observational 
techniques iYoachim & Dal canton 2006). Results from this simu- 
lation support the notion that thick disks are fo rmed during majo r 
mergers and the early assembly of galaxy disks jBrook et alj2004h . 
Furthermore, the last major merger should leave a clear signature 
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in the age distribution of thick disk (which will be older than the 
merger) and thin disk stars (mostly formed after the merger). 

In this realization no significant component of the stellar disk 
formed by accretion of stars in smaller satellites through minor 
mergers. While a substantial galaxy to galaxy scatter is expected, 
this is consist ent with results from our larger sample of simulated 
disk galaxies ( iBrooks & et~aill2008h and, compared with the exist- 
ing literature, a consequence of the smaller stellar masses of the 
gala xy satellites resu lting from the realistic feedback implementa- 
tion jGovemato et al. 2007). We verified that the bulge grows only 
modestly in mass from z=0.8 to z=0 and that no strong bar insta- 
bilities form after the last major merger Similarly the halo compo- 
nent shows minimal mass growth during the same time period. A 
close examination of the time evolution of the system shows sev- 
eral passages and disruption of small dark satellites, but they do not 
prevent the system from reforming a thin stellar disk. 



6 CONCLUSIONS 

We have analyzed an SPH, fully cosmological simulation of the 
evolution of a galaxy that grows an extended thin stellar disk after 
a major (1.6:1) merger at z=0.8. The disk dominates the light dis- 
tribution 3.5 Gyrs after the merger. By the present time the galaxy 
shares several properties with the observed population of disk dom- 
inated L* galaxies. 

This result is particularly relevant as the assembly history of 
the galaxy's parent DM halo is different from many previously pub- 
lished studies of cosmological simulations of disk galaxies, as it un- 
dergoes a late major merger, whereas previous works tended to se- 
lect galaxies with a relatively quiet merger history. The two progen- 
itors underwent major mergers themselves, at z ~ 2. The end result 
of this simulation strongly contradicts the notion that disk forma- 
tion requires a "quiet" halo merging history. In fact, the merging 
history of the halo picked for this study has often been considered 
hostile to the formation of extended disks at the present time. Our 
study provides strong support to the notion that galaxy disks can 
reform in a few Gyrs after a gas rich major merger, while a fraction 
of the pre-existing stellar disk can survive, even if faded by age and 
thickened by the strong interaction. 

The results of our analysis are made particularly robust by 
measuring the properties of the remnant light distribution rather 
than just that of the stellar mass. This approach is crucial, as it high- 
lights the effects of fading of the light from older stellar popula- 
tions and plays a major role in quantifying the predominance of the 
newly formed (and hence bluer and brighter) stellar disk. The main 
structural parameters of the galaxy at z=0 suggest that it has an Sa 
or Sb morphology: reddened i band B/D ~ 0.65, Rd = 7.2 kpc, g — r 
= 0.5, a W20/2 HI velocity width of 238 km/sec. Moreover, being 
able to measure the reddened colors and brightness evolution of the 
progenitors allows us a comparison with high-z galaxies, showing 
that they are representative of the population of blue, gas rich and 
moderately star forming galaxies observed at z>l. Verifying that 
the progenitors have some of the observed properties of high red- 
shift galaxies is important, as it makes the present day properties of 
the merger remnant more relevant, having been built from realistic 
progenitors. 

As our simulation includes a full treatment of the cosmolog- 
ical environment it includes a realistic treatment of a number of 
hydrodynamical processes that are necessary for the regrowth of 
stellar disks and that have to be simultaneously included: continu- 
ous inflow of gas from the cosmic web, stellar feedback, and cool- 



ing from hot halo gas. This explains the difference between our 
results and those of the collisionless simulations of IPurcell et al.l 
(2008), which showed significant disk heating due to interactions 
and accretion events. Lack of gas infall prevented the disks in their 
simulations to regenerate a new thin component (gas resupply was 
in fact advocated as a possible solution). In our simulation the ef- 
fect of interactions with infalling satellites, both luminous and dark, 
is naturally included. However, interactions and minor mergers are 
not strong enough to significantly disrupt or thicken the disk as it 
reforms from new gaseous material. 



Idealized hy drodyna mic simulations of b inary galaxy merg- 
ers teobertson e t al. 2006; Hopkins et al.ll2008l) did not include the 
continued infall of gas from cold flows and the hot halo. These 
works pointed out that without any subsequent accretion/growth 
onto the disk after the merger, a cold gas fraction in excess of 50% 
would be required in the disks of the progenitors in order to re-build 
a disk dominated galaxy after the merger. Here we have shown that 
such high disk gas fractions are not necessary, as cold flows and 
cooling from a hot halo make disk regrowth possible even at low 
redshifts, when the gas fraction in the progenitors' disks and the 
remnant is only ~20% at any given time. The analysis of the build 
up of the disk in our simulated galaxy highlights the dominant role 
of cold flows at high redshift. Cold flows funnel gas to the central 
disk on a time scale a few Gyrs shorter than gas that is first shocked 
to the virial temperature of the host halo and then cools onto the 
disk, leading not only to early disk star formation, but the creation 
of a large reservoir of cold gas. Thus, if a substantial fraction of this 
cold gas reservoir survives a merging event, it can provide a faster 
accretion rate onto the reforming disk than that available from cool- 
ing the hot gas in galaxy halos alone. Rapid disk reformation will 
also be aided if cold gas accretion onto the galaxy halo is still oc- 
curring, while gas cooling from the hot phase forms just 30% of the 
post merger stellar disk. However, this young component is more 
evident when the system is observed in bluer bands. It is important 
that these processes of disk formation and destruction be carefully 
implemented in analytical models that study the properties of large 
sample of galaxies. 



More work is also needed to make detailed quantitative pre- 
dictions about the morphology of galaxies formed in cosmologi- 
cal simulations and to make the results of our study more general. 
The maximum halo mass and lowest redshift at which mergers can 
regrow disks are obviously a function of th e feedback efficiency 
dBrook et alj|2004l : IScannapieco et alj|2008 l) and could therefore 
provide useful qualitative tests of models of SF and feedback. In 
dense environments, such as groups or clusters, the gas reservoir 
associated with cold flows and cooling halo gas will likely be dis- 
rupted by tidal forces and ram pressure stripping. Hence the disk 
re-growth process should be much less efficient, as expected by 
the observed correlation between galaxy morphology and environ- 
ment density. Also, higher resolution cosmological simulations will 
have to address the role of secular processes on the detailed struc- 
ture of the bulge and thin di sk components and their role in setting 
the B/D ratio of galaxies ( Debattist a et alj|2004l : iGenzel & et al I 
2008; Weinzirl et al. 2008). Still, results of the work presented here 
greatly alleviate the problem posed for the existence of disk galax- 
ies by the high likelihood of interactions and mergers for galaxy 
sized halos at relatively low z. 
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